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Abstract 

Recent measurements of jet cross sections in neutral-current and charged-current deep 
inelastic ep scattering are presented. The results of the QCD analyses on these measure¬ 
ments to determine the strong coupling constant are also reported. 


1 Introduction 

Jet production in neutral-current (NC) and charged-current (CC) deep inelastic ep scattering 
(DIS) provides a test of perturbative QCD (pQCD) calculations and of the electroweak sector of 
the Standard Model (SM). Jet cross sections allow the determination of one of the fundamental 
parameters of QCD, the strong coupling constant as, and help to constrain the parton densities 
in the proton. New particles or interactions may be observed by deviations of the measured jet 
cross sections with respect to the predictions. 

Up to leading order (LO) in as, jet production in NC and CC DIS proceeds via the quark- 
parton model (QPM) (Vq ^ q, where V = 'j, or W^), boson-gluon fusion (BGF) {Vg qq) 
and QCD-Compton (QCDC) (Vq ^ qg) processes (see hgure |^). 


eVv eVv eVv 



Figure 1: Quark-parton model (a), boson-gluon fusion (b) and QCD Compton (c) diagrams. 

The jet production cross section is given in pQCD by the convolution of the parton densities 
in the proton and the subprocess cross section, 


= X] dx fa{x,fip) daa{x,as{fiR),fi\,fiF), 
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where x is the fraction of the proton momentum taken by the interacting parton, fa are the 
parton distribution functions (PDFs) in the proton, fip is the factorisation scale, da is the 
subprocess cross section and fip is the renormalisation scale. 

The proton PDFs are determined from global hts to experimental data (see e.g. ]I|, H) and 
parametrise the long-distance structure of the target hadron. The subprocess cross section is 
calculable in pQCD at any order and represents the short-distance structure of the interaction. 

The determination of as from measurements of cross sections in DIS depends on the proton 
PDFs. Thus, to perform a QCD analysis, a functional form of the PDFs must be assumed 
and the determination of as should be made from observables with a small dependence on 
the PDFs, e.g. ratios of observables. For precise measurements of as, the experimental and 
theoretical uncertainties should be small. Small experimental uncertainties are obtained by 
measuring e.g. ratios of observables and the theoretical uncertainties are reduced by including 
next-to-leading order (NLO) QCD corrections to the cross section calculations. 


2 NLO QCD calculations 

Several programs are available to make NLO QCD calculations of jet cross sections in DIS: 
DISENT [|], MEPJET §, DISASTER++ § and NLO JET §. The NLO corrections include 
virtual corrections with internal particle loops and real corrections with a third parton in the 
hnal state. The programs differ in the treatment of the real corrections: DISENT, DISAS- 
TER-I--I- and NLOJET use the subtraction method whereas MEPJET uses the phase space 
slicing method. DISENT and DISASTER-I-+ are found to agree at the 2% level, and NLOJET 
is in good agreement with them. MEPJET agrees well with the other programs at LO, but 
presents differences of the order of 5% at NLO. MEPJET is the only program that includes 
and W exchange and only DISASTERS--!- allows the possibility to turn on the number of 
active flavours as a function of the scale. 

The calculation of jet cross sections at NLO in DIS depends on two scales, and fip. Two 
possible choices for these scales that have been considered in the analyses presented here are 
Q = v^, where is the virtuality of the exchanged boson, and the jet transverse energy. 

Since the NLO QCD calculations are for jets of partons and the measurements are done at 
the hadron level, the calculations need to be corrected for hadronisation effects. These effects 
have been estimated using Monte Carlo models for parton radiation and hadronisation |^, . 

The uncertainties of the QCD calculations include that due to terms beyond NLO, which is 
usually estimated by varying the renormalisation scale by factors between 1/2 and 2; it amounts 
to 5 — 10% and translates into ~ 8% in the determination of as{Mz)- The uncertainties on 
the value of as{Mz) and of the proton PDFs amount to ~ 5%. The size and uncertainty of 
the hadronisation corrections have been taken into account: the hadronisation correction factor 
amounts to ~ 5% with an uncertainty of the order of 2%. 


3 Inclusive jet cross sections 

The use of inclusive jet cross sections in a QCD analysis presents several advantages: inclusive 
jet cross sections are infrared insensitive and better suited to test resummed calculations and 
the theoretical uncertainties are smaller than for dijet cross sections. 

Inclusive jet cross sections in NC interactions have been measured in the Breit frame using 
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the kx cluster algorithm |^| in the longitudinally invariant inclusive mode for the kinematic 
region of > 125 GeV^, where is the virtuality of the exchanged boson. These cross 
sections refer to jets of transverse energy measured in the Breit frame, 35 above 8 GeV and 


jet pseudorapidity also in the Breit frame in the range —2 < < 1.8. The inclusive jet cross 

sections as a function of and i?^g are shown in hgure 
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Figure 2: Inclusive jet cross sections in NG DIS as a function of (a) and i^^g (c). Ratio 
data/theory as a function of (b) and (d). 


The measured cross section has a steep fall-off of hve (four) orders of magnitude over the 
measured range of (£';j^g). NLO QGD calculations using DISENT, the MRST99 0 set of 
proton PDFs and different choices of have been compared to the data. The calculations 
describe reasonably well the (-^tb) dependence of the cross section for > 500 GeV^ 
(E^'g > 15 GeV). 


Inclusive jet cross sections have been measured in the laboratory frame in NG [11 and GG 


interactions to search for deviations from the SM. The jets have been reconstructed using 


the kx cluster algorithm in the kinematic region of > 125 (200) GeV^ and at least one jet 
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-1 < < 2 for NC (CC) interactions is required. 


of > 14 (8) GeV and 

The inclusive jet cross sections as a function of for NC and CC interactions are shown 
in hgure ^ The cross section in NC interactions displays a steep fall-off of four orders of mag¬ 
nitude over the measured range. The behaviour of the distribution in CC interactions is 
very different: it is approximately constant at low E^^^ and falls less rapidly than in NC, and 
it approaches the NC cross section for E^^^ ~ 80 GeV. This behaviour can be interpreted as 
due to the presence of a massive propagator in CC events, as it has already been observed in 
measurements of the dependence of the inclusive CC DIS cross section |T^. The cross sec¬ 
tions for jets give an independent measurement of the different electroweak boson propagators 
in NC and CC DIS. 

The predictions of Monte Carlo models using either the color-dipole model (CDM) |0] or 
hrst-order QCD matrix elements plus parton-showers (MBPS) ||^ have been compared to the 
measurements. The calculations describe reasonably well the measured cross sections. No 
signihcant deviation from the SM is observed in NC or CC interactions up to the highest E^^^ 
measured value (~ 100 GeV). 


ZEUS 



Figure 3: (a) Inclusive jet cross section in NC and CC DIS as a function of E^^^. (b) Ratio 

data/theory as a function of in CC DIS. 


4 Dijet cross sections 


Differential cross sections for dijet production in NC DIS have been measured in the Breit 
frame for 470 < < 20000 GeV^. In this kinematic region, the experimental uncertainties 

on the reconstruction of both the positron and the hadronic hnal state are smaller than at 
lower Q^. In addition, the theoretical uncertainties due to the modelling of the hadronic hnal 


4 









state, to the proton PDFs and to the higher-order contributions are minimised. For these 
measurements, the kx cluster algorithm was run in the longitudinally invariant inclusive mode 
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the Breit frame. Two jets of > 8 GeV, > 5 GeV and —1 < < 2 in each 

event were required, where is the transverse energy of the jet in the Breit frame 

with the highest (second highest) transverse energy in the event. The use of asymmetric cuts 
on B avoids infrared sensitive regions where the behaviour of the cross section as predicted 
by the NLO QGD programs is unphysical. 

The differential dijet cross sections as functions of where jet 1 (2) is the jet with 

highest (lowest) pseudorapidity in the Breit frame, and are presented in hgure |[ NLO QGD 
predictions calculated using DISENT are compared to the measurements. The predictions, 
which assume as{Mz) = 0.118, provide a good overall description of both the shape and 
magnitude of the measured cross sections. The dijet fraction, R 2 + 1 , dehned as the ratio of the 
dijet to the inclusive cross section (see hgure §(a)), has been measured as a function of Q^. The 
dijet fraction increases with increasing due to phase-space effects. For the cross sections as 
a function of the jet transverse energies and Q^, there is agreement at the ~ 10% level between 
data and theory over four orders of magnitude, demonstrating the validity of the description of 
the dynamics of dijet production by the NLO QGD hard processes. 

Dijet cross sections in GG interactions have been measured [|T^ in the laboratory frame 
as a function of the transverse energy of each of the two highest jets (hgure |^(a)) and 
as a function of the invariant mass, mjj, of these two jets (hgure §|(b)). The jets have been 
reconstructed using the kx cluster algorithm in the kinematic region of > 200 GeV^ and at 
least two jets of F;)?* > 8 GeV and — 1 < < 2 are required. Values as high as 60 GeV for 

E^r^^ and mjj are accessible with the luminosity used in this analysis. The shapes of the F;)^* 
distributions at low F;)?* in inclusive jet (see hgure ^ and dijet production are very diherent: 
this is interpreted as the inclusive jet cross section being governed by electroweak interactions 
and the heavy mass of the boson propagator at low F^*, whereas the dijet cross section is 
mainly dictated by QGD. The predictions of the GDM and MEPS models describe reasonably 
well the shape of the measured dijet cross sections, but they underestimate the normalisation 
by 50%; this is understood since these models only include up to 0{as) diagrams. 

5 Jet substructure 


The investigation of the internal structure of jets gives insight into the transition between a 
parton produced in a hard process and the experimentally observed spray of hadrons. At suf- 
hciently high F;)?*, where fragmentation effects become negligible, the jet structure is expected 
to be calculable in pQGD. The lowest non-trivial order contribution to the measurements of jet 
substructure is given by 0{aas) calculations. Thus, measurements of jet substructure provide 
a stringent test of pQGD and allow a determination of as by comparing NLO calculations to 
the measurements. At present, this is only possible for jets dehned in the laboratory frame 
since only in this frame can three partons be inside one jet (see hgure 

The jet substructure was studied by means of the integrated jet shape in inclusive jet pro¬ 
duction in NG interactions [O using the kx cluster algorithm in the laboratory frame. The 


integrated jet shape is dehned as the average fraction of the jet’s transverse energy that lies 
inside a cone in the rj — ip plane of radius r concentric with the jet axis. 
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Figure 4: Dijet cross sections in NC DIS as functions of (a) (b) cind (c) . In 

(c), the inclusive differential cross section (datot) is also shown as a function of . 


where ET{r) is the transverse energy associated with the jet within the given cone of radius r 
and iVjets is the total number of jets in the sample. Only the particles assigned to the jet are 
considered. The integrated jet shape as a function of the radius r in different regions of is 
shown in hgure p. 

In pQCD, (1 —'0(r)) is calculated as the fraction of the jet’s transverse energy, due to parton 
emission, that lies outside of the cone of radius r. 


(1 -^(r)) = 


f (IEt Et [da{ep 2 partons)/di?T] 


4 "* a. 


jet 


( 4 ") 


where crjet(F'y*) is the cross section for inclusive jet production. NLO QCD predictions for 
the integrated jet shape are derived from the above formula by computing the numerator to 
O(aag) and the denominator to O(aas). Calculations using DISENT have been compared to 
the data. They give a very good description of the measurements for r > 0.2 (see hgure H). 

The jet substructure was also studied by means of the mean subjet multiplicity in NC [jT^ and 
CC 0] interactions. Subjets are jet-like objects within a jet which are resolved by reapplying 
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Figure 5: (a) The dijet fraction, R 2 + 1 , in NC DIS as a function of Q^. (b) The as values 

determined from the QCD fit of the measured dijet fraction R 2+1 as a function of Q. 


the kx cluster algorithm until for every pair of particles the quantity 


dij = min(ET,j, - VjY + 

is above j/cut ■ where Dent is the resolution parameter. 

The mean subjet multiplicity as a function of the resolution parameter Uent is presented in 
figure I for NC interactions. In pQCD, (usubjet — 1) is calculated as the ratio of the cross section 
for Usubjet — 1 subjets over that of inclusive jet production, 


('^subjet 1) — 


cr, 


Flsubjet 


-i(4“) 




NLO QCD predictions for the mean subjet multiplicity are derived from the above formula by 
adding 1 and computing the numerator (denominator) to 0{aag) {0{aas)). NLO calculations 
computed using DISENT and different parametrisations of the proton PDFs, which were de¬ 
termined assuming different values of as, have been compared to the data. The NLO QCD 
calculation using the CTEQ4M proton PDFs and as{Mz) = 0.116 gives a good description of 
the measurements. 

The measured mean subjet multiplicity as a function of i/cnt for different regions in and 
is shown in hgures and respectively, for CC interactions. The predicted mean subjet 
multiplicities at the hadron level using the CDM and MEPS models have been compared to the 
measurements. The measured mean subjet multiplicities are well described by the predictions. 

The 4* dependence of the integrated jet shape at a hxed value of r = 0.5 (hgure |^(a)) and 
that of the mean subjet multiplicity at a hxed value of i/cut = 10“^ (hgure |^(b)) in NC interac¬ 
tions have also been studied |^. The jet shape increases and the subjet multiplicity decreases 
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Figure 6: Dijet cross sections in CC DIS as a function of (a) and mjj (c). Ratio data/theory 
as a function of (b) and mjj (d). 




Figure 7: Diagrams in O(aa^) calculations. 


as 77;^* increases: the jets become narrower as increases. NLO QCD calculations using 
DISENT and different parametrisations of the PDFs have been compared to the data. They 
give a good description of the measurements and display the sensitivity of these observables to 


as- 

Figure |^(a) shows the measured (rigubjet) at 2/cut = 10“^ as a function of in CC inter¬ 
actions. It exhibits no signihcant dependence on as in the case of the jet shape [^. The 
measured (usubjet) at pcut = 10“^ as a function of E^*' is presented in hgure |T^(b). It decreases 
as E^r^^ increases and thus the jets become narrower as E^^ increases. In both hgures, the 
predictions from the CDM model are in good agreement with the data whereas the predictions 
from the MEPS model are slightly above the data. In hgure p!B|(b), the measurements from 
hgure |T3(b) are also shown. The measured (usubjet) at ycut = 10“^ in CC and NC are very 
similar. This similarity can be attributed to a large content of hnal-state quark jets in these 
two processes and shows that the pattern of QCD radiation close to a primary quark is to a 
large extent independent of the hard scattering process. 
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Figure 8: Integrated jet shape as a function of r in different regions of in NC DIS. 


6 QCD analyses 


6.1 Determination of as 

In recent publications |^, |^, a new method to determine as has been presented. The 
method consists of performing NLO calculations of an observable using the hve CTEQ4 sets 
of the “A-series” |^; the value of as{Mz) in each calculation is that of the corresponding set 
of PDFs. The calculations for each as{Mz) are used to parametrise the as dependence of the 
observables, and then, from the measurements, a value of as{Mz) is directly extrapolated. The 
measurements of the inclusive jet cross section P] and the dijet fraction |]T^ as a function of 
and the integrated jet shape [|T^] and mean subjet multiplicity as a function of E', 
been used to determine a value of as{Mz 


have 


From each measurement in a given Q‘ 


or 


region, a value of as{Mz) was determined. 


The results are shown in hgure separately for the inclusive jet cross section, the mean subjet 
multiplicity and the integrated jet shape. 

A combined value of as{Mz) from each observable has been obtained by performing a 
£t to the data in the following kinematic ranges: > 500 GeV^ for the inclusive jet cross 

section, > 470 GeV^ for the dijet fraction, > 21 GeV for the integrated jet shape and 
> 25 GeV for the mean subjet multiplicity. These are the kinematic regions where the 
hadronisation correction factors to the NLO calculations are small for each observable. The 
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combined values of as{Mz) determined from each observable are: 


inclusive jet cross section: 
dijet fraction: 
integrated jet shape: 
mean subjet multiplicity: 


as{Mz) = 0.1190 ±0.0017 (stat) toS (exp) (th); 

as{Mz) = 0.1166 ±0.0019 (stat) (exp) (th); 

as{Mz) = 0.1179 ±0.0014 (stat) (exp) (th); 

as{Mz) = 0.1185 ±0.0016 (stat) IqS (exp) +o;oo 7 ? (th). 


The values of as{Mz) obtained from the different observables 
and with the world average [^, [T^ (see hgure |T^). 


are compatible with each other 


6.2 Energy scale dependence of as 

The QCD £t of the dijet fraction has been repeated in hve bins to test the scale dependence 
of the renormalised strong coupling constant. The method of the £t is the same as outlined 
above, but the as dependence of the dijet fraction was parametrised in terms of as{{Q)), where 
{Q) is the mean value of Q in each bin. The measured as{{Q)) values are shown in figure ^(b). 
The measurements are compared with the renormalisation group predictions obtained from the 
PDG as{Mz) value and its associated uncertainty. The values are in good agreement with the 
predicted running of the strong coupling constant over a large range in Q. 
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Figure 10: Mean subjet multiplicity as a function of t/cut in different regions in CC DIS. 

7 Summary and conclusions 

Recent measurements of jet cross sections in neutral-current and charged-current deep inelastic 
ep scattering have been presented. From QCD analyses of these measurements, determinations 
of the strong coupling constant as have been obtained. The results are consistent with the 
world average value of as{Mz)- The measurements have been performed in regions of phase 
space where the experimental uncertainties are small and thus the determinations are precise. 

However, the next-to-leading order (NLO) calculations used in the QCD analyses have rela¬ 
tively large uncertainties. For more accurate measurements of the parameters of QCD, improved 
calculations are needed. For instance, higher-order corrections are required to reduce the largest 
theoretical uncertainty, namely that of the terms beyond NLO. 

Acknowledgments. This report is partially based on a talk given at the International Eu¬ 
rophysics Conference on High Energy Physics, Budapest, Hungary, July 2001. I would like to 
thank J. Terron for a critical reading of the manuscript. 


References 

[1] H.L. Lai et al, Phys. Rev. D 55 (1997) 1280. 

[2] A.D. Martin et al, Eur. Phys. Jour. C 14 (2000) 133. 

[3] S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291; Erratum B 510 (1997) 503. 


11 




ZEUS 



^cut 


^cut 


Figure 11: 


Mean subjet multiplicity as a function of ycnt in different regions in CC DIS. 


[4] E. Mirkes and D. Zeppenfeld, Phys. Lett. B 380 (1996) 205; E. Mirkes, TTP-97-39 (|hep- 
I ph/9711^ . 


[5] D. Graudenz, Ihep-ph/ 9710244 . 

[ 6 ] Z. Nagy and Z. Trocsanyi, Piep-ph/0104315| . 

[7] Y. Azimov et al, Phys. Lett. B 165 (1985) 147; G. Gustafson, Phys. Lett. B 175 (1986) 
453; G. Gustafson and U. Petersson, Nucl. Phys. B 306 (1988) 746; B. Andersson, G. 
Gustafson and L. Lonnblad, Z. Phys. C 43 (1989) 625. 


[ 8 ] G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm. 101 (1997) 108. 

[9] ZEUS Gollaboration, “Inclusive jet cross sections in neutral current deep inelastic scatter¬ 
ing in the Breit frame and determination of 0:5 at HERA ”, paper 637, submitted to the 
International Europhysics Gonference on High Energy Physics 2001, Budapest (2001). 


[10] S. Gatani et al, Nucl. Phys. B 406 (1993) 187; S.D. Ellis and D.E. Soper, Phys. Rev. D 48 
(1993) 3160. 


[11] ZEUS Gollaboration, “Inclusive jet cross sections in deep inelastic ep scattering at high 
paper 536, submitted to the International Europhysics Gonference on High Energy 
Physics 1999, Tampere (1999). 


[12] ZEUS Gollaboration, “Measurements of jet cross sections and jet substructure in charged 
current deep inelastic scattering at HERA”, paper 646, submitted to the International 
Europhysics Gonference on High Energy Physics 2001, Budapest (2001). 


12 















<+(r=0.5)> 



Figure 12: Integrated jet shape at a fixed value of r = 0.5 (a) and mean subjet multiplicity at 
a fixed value of ?/cut = 10“^ (b) as a function of in NC DIS. 


[13] ZEUS Collaboration, J. Breitweg et al, Eur. Phys. Jour. C 12 (2000) 411. 

[14] ZEUS Collaboration, J. Breitweg et al, Phys. Lett. B 507 (2001) 70. 

[15] ZEUS Collaboration, “Measurements of jet substructure in neutral current deep inelastic 
scattering and determination of as at HERA”, paper 641, submitted to the International 
Europhysics Conference on High Energy Physics 2001, Budapest (2001). 

[16] ZEUS Collaboration, J. Breitweg et al, Eur. Phys. Jour. C 8 (1999) 367. 

[17] Particle Data Group, D.E. Groom et al, Eur. Phys. Jour. C 15 (2000) 1. 

[18] S. Bethke, J. Phys. G 26 (2000) R27. 


13 










sbj 


ZEUS 


ZEUS 




Figure 13: Mean subjet multiplicity at a fixed value of i/cut = 10 ^ as a function of (a) 
and (b) in CC DIS. 


ZEUS ZEUS 




a. (Mz) 


a, (Mz) 



world 

average 

<- 


30 40 50 60 70 


Ef [GeV] 


Figure 14: Extracted values of as{Mz) from the inclusive jet cross section (a), the mean 
subjet multiplicity (b) and the integrated jet shape (c) in NC DIS. 
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Figure 15: A compilation of recent determinations of as{Mz) using jet production in NC DIS. 
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